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ABSTRACT: We calculate the free energy and the pressure of a weakly slip- 
linked Gaussian polymer chains. We show that the equilibrium statistics of a 
slip-linked system is different from one of the corresponding ideal chain system 
without any constraints by slip-links. It is shown that the pressure of a slip- 
linked system decreases compared with the ideal system, which implies that 
slip-linked chains spontaneously form aggregated cluster like compact struc- 
tures. These are qualitatively consistent with previous theoretical analyses 
or multi chain simulations. We also show that repulsive potentials between 
chains, which have been phenomenologically utilized in simulations, can cancel 
the artificial pressure decrease. 
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INTRODUCTION 

The motion of polymer chains in melts or solutions are strongly constrained due to the 
entanglement effect \ if the degrees of polymerization are large. Although the entanglement 
affects dynamic properties drastically, it does not affect the static properties. For example, 
in polymer melts, conformations of polymer chains are well described by the ideal, Gaussian 
statistics^ beyond a certain length scale. Therefore, if we are interested only on static 
properties, the Gaussian chain model can be reasonably utilized. Actually theories based on 
the Gaussian statistics achieved great success to predict static phase behaviours of polymer 
blends or block copolymers However, once we consider dynamics of entangled chains, the 
situation becomes complicated so much. 

It is informative to consider simple toy models which show non-trivial dynamical be- 
haviours due to topological constraints. It was reported that molecules composed of infinitely 
thin rods exhibit glass transition by increasing the density-"-, while thermodynamically they 
are ideal gases. The situation is almost the same for entangled polymers. The entanglement 
effects exist even for ideal chains as long as they cannot cross each other. (Although in- 
teractions between monomers or between monomer and solvent alter rheological behaviours 
somehow, such effects are secondary.) Unfortunately, it is not easy to treat such dynamic 
effects theoretically. Thus several phenomenological models such as the tube model- or 
the slip-link model- have been proposed and utilized. The slip-link model represent an 
entanglement as a linking point on polymers which can slide along the chains (slip-links). 
So far, various slip-link type models^"— have been utilized to study rheological properties 
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of entangled polymers, and slip-link models have reproduced various rheological properties 
well. 

However, it is reported that the slip-link models can exhibit unintuitive statistics in 
several situations^""—. The simplest example is the statistics of a single Gaussian chain with 
slip-links. Rieger— ~— calculated the partition function of a Gaussian chain with a slip-link 
and showed that this system can exhibit spontaneous symmetry breaking and bifurcation. 
Metzler, Kardar and co workers^"— studied a single ring polymer with slip-links as a model 
of a knotted ring polymer, and found that slip-linked portion localizes on a ring. This is 
similar to Rieger's results. They reported that the localization behaviour is observed even 
in presence of the excluded volume interaction. Sommer^^ studied two chains bound by a 
slip-link, and found similar unintuitive behaviours. He also found that when multiple chains 
are bound by multiple slip-links, slip-links can form compact aggregated structures. These 
results imply that various slip-link models may exhibit unphysical results. 

We consider that macroscopic statistical properties of slip-linked chains may be affected 
by slip-links. This can be serious problems when we utilize slip-link models to study var- 
ious macroscopic properties (such as pressure or rheology). As far as the authors know, 
so far studies for such unintuitive statistics are limited for simple systems such as single 
chain systems. How slip-links affect macroscopic thermodynamic properties is not clear. In 
principle, if the statistics of the system is explicitly specified (if the equilibrium probability 
distribution is explicitly given), we can check whether the macroscopic thermodynamic prop- 
erties are really affected or not. For example, in single chain slip-link models with the mean 
filed type approximation— or in virtually coupled many chain slip-link models^, we can 
straightforwardly show that there are no unintuitive behaviors (as discussed in Appendix 
A). 

The situation is not simple when we consider multi chain slip-link models^. In this 
work we consider statistics of weakly slip-linked multi polymer chains theoretically. We 
calculate the partition function of weakly slip-linked Gaussian chains by utilizing the grand 
canonical ensemble type treatment developed by Schieber^fl. We show that we can calculate 
the explicit forms of partition function and free energy as perturbation expansions. The 
free energy of weakly slip-linked chains are shown to be different from one for ideal chains 
without any constraints by slip-links. We also show that slip-links decrease the pressure, 
and thus repulsive interactions between chains are required for such systems to recover (or 
mimic) the ideal statistics. 

THEORY 

Thermodynamic Equilibrium Distribution of Gaussian Chains 

Before we calculate the partition function for slip-linked chains, we first consider the partition 
function of Gaussian chains without any constraints by slip-links. (We call such systems as 
"the ideal chain systems" or "the ideal systems", in the foUowings.) As we mentioned, 
even if chains are well entangled, equilibrium statistical properties can be well described 
by the Gaussian chain model since the entanglement effect is purely kinetic. As pointed 
simply by Doi^^, the entanglement effects exist even for ideal chains without excluded volume 
effects. Therefore we limit ourselves to ideal chains and ignore interaction potentials between 
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monomers although they may not be neghgible for some cases (especially for small length 
scales). Thus the partition function or the free energy for ideal chains can be directly 
applied both to unentangled or entangled polymers. Although one may consider this is an 
oversimplification, as we will show in the foUowings, we can still obtain nontrivial results 
for slip-linked chains. We will compare the results of this subsection with the statistics of 
slip-linked chains later. 

Some expressions in this subsection can be utilized in the next subsection, to calculate 
the partition function for slip-linked chains. Thus here we calculate the partition function 
for Gaussian chains rather verbosely. It is convenient to start from a single chain partition 
function, since the partition function of the system can be easily calculated from the single 
chain partition function in non-interacting multi chain systems. We describe the conforma- 
tion of a polymer chain as R{s) (0 < s < A^). s is the index along the chain and N is the 
degree of polymerization. The single chain partition function without any constraints can be 
formally described by using the Edwards Hamiltonian and the functional integral expression. 
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where ?{o[-R(')] is the Edwards Hamiltonian and J VR means the functional integral over 
-R. We assume that the functional integral j VR contains an appropriate measure so that 
the resulting functional integral becomes dimensionless. The subscript 1 in the left hand 
side of eq ([T]) indicates that it is the partition function for a single chain. The Edwards 
Hamiltonian for a Gaussian chain is expressed as 
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Here ks is the Boltzmann constant, T is the temperature, and b is the segment size. (As we 
mentioned, the interaction between monomers are not considered in this work and thus we 
have no interaction energy in eq ([2]).) 

We introduce the constraint which fixes the position of the chain at the index si, -R(si), 
to ri. From translational symmetry (translational invariance), the partition function under 
this constraint is independent of Si and Vi. Then we can define the following partition 
function. 



Q= VR [A=^5(ri - R{s,))] exp 



ksT 



(3) 



Here we introduced the dimensional factor A'^ (with A being the thermal de Broglie wave 
length of the chain) to make the partition function dimensionless. Eq ([3]) corresponds to 
the partition function for the internal degrees of freedom. By using eq (I3l), the single chain 
partition function ([1]) can be rewritten as 
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where V is the system volume. 
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Next we consider a canonical ensemble of Gaussian chains. The partition function for M 
Gaussian chains can be simply written as follows, by using the single chain partition function 

M! M!A3A^ 

The superscript (0) is used to indicate an ideal system. The free energy of the system is 
simply calculated to be 



jr(o) = -keTXTiZ^^^ = keTM 
Thus the pressure becomes as follows. 
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In -TTT^ - 1 
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(6) 



(0) ^ ^ ksTM 

dV V ^ ^ 

This is nothing but the standard van't Hoff form (the equation of state for ideal gas). The 
pressure is simply proportional to the polymer chain density M/V. We again emphasize 
that eq ([7]) holds even if chains are entangled, because the entanglement effect is purely 
kinetic. Although results shown here are rather trivial, it is not trivial whether the slip- 
linked chains have the same statistics or not. In the following sections, we calculate the 
static thermodynamic properties of slip-linked chains. 



Partition Function of Slip-Linked Chain Clusters 

We consider weakly slip-linked chains in this and following subsections. For simplicity, we 
assume that slip-links do not interact each other, and a single slip-link always bind two chain 
portions in space. Under the constraints by slip-links, the thermodynamic properties of the 
system will be altered. It can be shown by calculating the partition function for a slip-linked 
chains. If we assume that the constraints by slip-links are sufficiently weak (or chains are 
sufficiently short), it is justified to neglect complicated network like structures with many 
slip-links. We show an image of a weakly slip-linked multi chain system in Figure [1] Then 
we need to calculate only several simple slip-linked structures, to obtain the expression of 
the partition function of the system. This is similar to the case of weakly interacting real 
gas^. If the interaction between gas molecules is sufficiently weak, we only need several 
lower order Mayer clusters to study thermodynamic properties. 

The most simple slip-linked structure is two chains paired by one slip-link. The schematic 
image of a chain-pair is depicted in Figure [2]^a). We may call such a structure as a "slip- 
linked chain cluster" or simply a "cluster", in analogy to the Mayer clusters. It is easy to 
show that the expression for the partition function is different from the ideal case. The 
partition function can be calculated straightforwardly by imposing the slip-linking effect by 
a delta function. We express the conformations of two chains by Ri{s) and R2{s). Two 
points on these chains, -Ri(si) and R2{s2) are coupled by a slip-link (si and S2 are indices 
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of slip-linked points). The partition function is calculated as follows. 
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where we introduced the subscript 2(1, 1) which indicates that the cluster is composed of 
two chains and each chain has one slip-linked point. 

If there are two slip-links, more complicated slip-linked structures can be formed. In this 
case, there are three different possible cluster structures (Figure [2]^b)-(d)). The cluster shown 
in Figure [2](b) is composed by three chains. The clusters shown in[2](c) and (d) are composed 
by two chains, but their topologies are different. After straightforward calculations, we have 
the following partition functions for these clusters. 
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The subscripts 3(1,1,2), 2(1,3) or 2(2,2) indicate the number of chains and slip-links on 
each chain, as before. For example. The cluster 3(1, 1,2) is composed by three chains and 
they have 1, 1, and 2 slip-linked points. 

There are other clusters which are composed of one or two slip-links. As discussed in 
Appendix B, these clusters have self-slip-linked structures and such self-slip-linked clusters 
have no effect on static thermodynamical properties. Intuitively, we can understand this as 
follows. In the self-slip-linked structures, slip-linked portions strongly localize^"—. Then the 
partition function for a self-slip-linked structure is almost the same as the partition function 
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for a structure without self-slip-links. Thus such clusters practically do not contribute to 
the partition function. Anyway, we do not need to consider them when we calculate the 
partition function of slip-linked chains. All the clusters considered in this subsection are non 
self-slip-linked clusters and contribute to thermodynamic properties. 



Equilibrium Free Energy for Weakly Slip-Linked Chains 



We can calculate the partition function of the system from partition functions for clusters 
fl8|)- f|TT]) . We start from a nearly ideal multi chain system where only small fraction of chains 
forms slip-linked clusters. Under such a condition, the system partition function can be 
calculated only by using the partition functions for a free chain and a cluster 2(1, 1). 

According to the Schieber's theory^S-, the distribution of slip-links can be regarded as 
a sort of ideal gas in a grand canonical ensemble, and they are controlled by the effective 
chemical potential for slip-links. Then the partition function of the system can be described 
as 
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where M is the total number of chains in the system and e is the effective chemical potential 
for a slip-link. The sum in eq ( IT2l) is taken for M\ and M2 under the constraint M\ -|- 2M2 = 
M. The calculation becomes easier if we consider a grand canonical ensemble of chains, 
instead of a canonical ensemble. If we introduce the chemical potential for a chain, /x, the 
grand partition function becomes 
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Eq ( !T3|l can be easily calculated. 

S = exp {e^l^^^Z^ + e^^'^+^^/'^^^Zsd,!)] 
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Thus we obtain the following expression for the grand potential J 

J=-kBT\n'E 

1 



KbI ^3 e 



1 + -e(M+f)/fcsT'g 



Now we can calculate the free energy by using the Legendre transform. 
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The pressure of the system is calculated to be 

kBTM 



P 



dV 



2V 



1 + 



2 V V 2 



-1 



Clearly the free energy fll7p and the pressure (1181) are different from ones for an ideal system 
(eqs © dZD) 

It is convenient to introduce the average number of slip-linked points on a chain, Z. The 
average number of slip-linked points (entanglements) on a chain can be calculated directly 
from the free energy. 

-1 



2 dJ^ 



1 



1 / ,^ ^MA^ 1 

2V V 2 



(19) 



Here the factor 2 is introduced because one slip-link corresponds to two entanglement points 
(see Figure [2]). From eq ( lT9l) we find that the average number of slip-links on a chain depends 
on the polymer chain density M/V. We can modify eq (fTOll and obtain the following relation. 
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Eq fl2Ul) gives the relation among the effective chemical potential, the chain number density, 
and the number of slip-links on a chain. Since we are considering weakly slip-linked chains, 
the number of slip-linked points on a chain should be sufficiently small {Z <^ 1) and ther- 
modynamic quantities can be expanded into power series of Z around Z = 0. The power 
series expressions will help us to interpret the physical meanings of the results. Eq (!20l) can 
be expanded as follows. 

^e/k,T^^^^ + 0{Z') (21) 

By using (12T|) . we can rewrite the free energy (fT7|) and the pressure (fT8|) simply as follows. 
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It is now clear that eqs fl22|) and (p3|) reduces to eqs ([6]) and ([7]), respectively, at the limit of 
Z — )■ 0. We can rewrite eqs (122|) and fl23l) as follows by using eqs ([6]) and ([7]). 



jr _ jr(o) = -htBTMZ + 0{Z^) 



P _ p(o) 



(24) 
(25) 



Eqs flMl) and fl25l) give the expressions for the excess free energy and the excess pressure. From 
eq fl2^ we find that the system free energy is decreased by slip-links. Eq (1251) means that the 
pressure of a weakly slip-linked chain system becomes slightly lower than the pressure of the 
corresponding ideal chain system. A weakly slip-linked system does not obey the equation 
of state for an ideal gas. 
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Effect of Higher Order Clusters 



In the previous subsection, we have considered only the first order cluster. In this subsection 
we consider the effect of the higher order clusters. If we retain the clusters up to the second 
order, the expression for the grand potential becomes 
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From the results of the previous subsection, we expect that ^ = e'^'^+'^^/'^^-^Q becomes a small 
perturbation parameter which represents the strength of constraints by slip-links. (^ may be 
interpreted as the effective fugacity for slip-links. ^ = corresponds to a non slip-linked, free 
ideal chain system.) Then thermodynamic quantities can be expanded into power series of 
^ around ^ = 0. Since all the clusters which contain two or less slip-links are included in eq 
( l26l) . we can calculate the expansion up to 0(^^). From eq ( l2Ti) we find that 0{Z) = O(^), 
and thus we can calculate the expansion series up to O(Z^). 
The number of chains in the system is calculated to be 
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Then the chemical potential can be related to M and C, as follows. 
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After the straightforward calculations, we have the following expressions for the grand par- 
tition function or the free energy. 



J = -keTM 



1 - -£ + -r 

2^ 3^ 



2(61-32^2) 2 ( ^ \ 



15 



3/2 



JT = jr{o) _ keTM 
The pressure can be obtained from eq (130|) . 

p = p(o) - keT 
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(30) 

(31) 



From eq (!3T|) it is clear that by increasing ^ (by increasing the strength of the slip-linking 
effect), the pressure decreases. 
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The average number of shp-hnks Z can be calculated in the same way as the previous 
case. 

Eq (l32i) means that the number of slip-links monotonically increases as ^ increases, but the 
dependence of Z to ^ is nonlinear. By inverting eq f p2|) . ^ can be expressed as a function of 
Z as follows. 



%~ 8(61 - 32v/2) ^, V- / 3 



(33) 



Finally the free energy ( l30l) and the pressure ( l3Ti) can be rewritten as follows, by using the 
average number of slip-links on a chain, Z . 
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Comparing eqs ([23]) and fl35p . we find that the pressure becomes slightly lower if we take 
account of the higher order clusters. Intuitively this looks inconsistent with eq fl35|) in which 
the pressure becomes slightly higher by the higher order clusters. This comes from the 
nonlinear relation between ^ and Z (eq (!33l)). From eqs ( l27|) and ( l32l) . Z can be expanded 
into power series of e*^/^^"^ and M/V. 
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+ 0(e^^/'=^^(M/V)^) 



Although eq ( l36l) is not simple, we can easily find that Z depends on the chain density M/V 
nonlinearly. 

DISCUSSION 



Here we discuss about the pressure of slip-linked multi chains in detail. As we have shown, 
the pressure of a slip-linked system decreases from the pressure of the corresponding ideal 
system. Because the thermodynamic control parameter for slip-linked systems is the chain 
density M/V and the effective chemical potential for slip-links e, we rewrite the pressure 
( 135|) to see how these parameters affect the pressure. From eqs ( l35i) and ( l36l) . we can rewrite 
the excess pressure as follows. 



„e/fcsT A 3 , 2(61 - 32y/2) 2e/kBT .6 f ^ \ 

^ 15 \27rNby 



3/2 -1 



y) +0{{M/Vf) 
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The leading order term of the excess pressure is proportional to (M/V^)^, and we can in- 
terpreter it as the second order virial. The second order virial in eq (!37|) gives a negative 
correction to the pressure. Roughly speaking, the negative second order virial gives the 
same contribution as attractive interaction between chains. Clearly there are no interaction 
in ideal systems, and thus the second order virial for an ideal system is exactly equal to zero. 
Therefore we can conclude that the negative second order virial is an artifact of the slip-link 
model. It should be noted that from eq fl37j) the excess pressure is proportional to fc^T. 
This means that the origin of the excess pressure is entropic. This can be interpreted that 
the entropy loss by slip-links cause the negative virial (or the effective attraction between 
chains) . 

This result implies that we should introduce an effective repulsive interaction between 
chains to cancel the artificial second order virial and recover the correct static thermody- 
namic properties. For example, we may introduce a phenomenological repulsive interaction 
potential between monomers. If we assume that the interaction is sufficiently weak, the 
resulting total pressure, P, can be formally expressed as 

P^P + ksTBA-^] (38) 



where B2 is the second order virial coefficient calculated from the interaction potential, and 
NM/V corresponds to the monomer density. By tuning the interaction potential, we can 
effectively cancel the contribution of the artificial second order virial. The potential should 
be tuned to satisfy the following condition. 



1 /,^^ 2(61-32v/2) ( 3 \ 

2 ^ 15 \2'KNhy 
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(39) 



Notice that, even if the second order virial is successfully cancelled, there are still higher 
order virials and the full ideal chain statistics cannot be reproduced in a slip-linked system. 
This situation is somehow similar to polymer chains in a ^ solvent or real gases at the Boyle 
point. Although the second order virial is cancelled in a ^ solvent, there still exists a non-zero 
third or higher order virials. Therefore generally it is impossible to fully restore the ideal 
chain statistics for a slip-linked system, by an artificial interaction potential. Similarly, we 
can introduce an effective interaction between slip-links. In this case, we have the following 
total pressure instead of eq (1551) . 



P^P + keTB^y-^j (40) 

Here B2 is the second order virial for the interaction between slip-links. By tuning the 
interaction potential, the total pressure can be reduced to the ideal pressure up to the 
second order in M/V , as before. 

Here it should be pointed that there is no such an artificial pressure decrease in the single 
chain type models, such as slip-linked single chain models by Schieber and coworkers^iii^^^, 
the dual slip-link model by Doi and Takimoto--, or the slip-spring model by Likhtman— . In 
these models, the partition function is essentially the same as the partition function of ideal 
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chains—. (See Appendix A.) It is reasonable to consider that the main reason why we have 
the artificial decrease of pressure in a multi chain slip-linked system is that slip-links bind 
two chains spatially. Chains bound by slip-links no longer follow the ideal statistics and thus 
the resulting thermodynamic properties become different from ideal ones. 

Let us consider a non-interacting slip-linked multi chain system again, from a different 
aspect. Analyses and discussions above are limited for weakly slip- linked systems. Here 
we consider strongly slip-linked systems by utilizing a simple approximation. Although in 
general it is impossible to obtain the analytical expressions for strongly slip-linked systems, 
if we consider only limited sets of slip-linked clusters (which have simple topologies), it 
is possible to obtain approximate expressions. One of the simplest approximations is to 
use only one dimensionally connected, chain-like clusters. By utilizing this chain clusters 
approximation, exact expressions for the free energy or the pressure can be obtained. Detail 
calculations are described in Appendix C. The pressure by the chain cluster approximation 
becomes as follows. 



where < Z < 2 is the average number of slip-linked points on a chain. Clearly, at the limit 



This can be interpreted that all the chains aggregate into a single slip-linked cluster, and 
thus this system is thermodynamically almost the same as a single ideal gas particle in a 
large box. The pressure of such a system is approximately zero, because the pressure is 
proportional to the number density of particles and the density is approximately zero. 

We expect that the homogeneous state becomes unstable and slip-linked chains sponta- 
neously form compact aggregates if the slip-linking effect becomes sufficiently strong. This 
is qualitatively in consistent with the theory of Sommer— . He predicted that slip-links ag- 
gregate even if there is no other direct interaction between chains. The compact aggregate 
formation in slip-linked systems is actually observed in multi chain slip-link simulations—. 
Masubuchi et al— introduced phenomenological osmotic pressure for their model to avoid 
aggregation. Several different forms of the phenomenological osmotic pressure have been 
utilized— >2^~"—, and it was found that all of these different forms give qualitatively the same 
result. So far, the use of such phenomenological osmotic pressure models are not theoretically 
justified. 

Our result indicates that osmotic pressure works as effective repulsion interaction between 
chains. If the osmotic pressure is sufficiently high, the artificial attraction by slip-links can 
be cancelled. Detail forms of the repulsive interaction potentials are not essential as long 
as they give appropriate positive contributions to the pressure. This is clear if we consider 
the condition for a weakly slip-linked system. There are various possible potentials which 
give the same second order virial and satisfy the condition (139|) . Therefore we conclude 
that the detail forms of osmotic pressure are not important in multi chain slip-link models. 
Our analysis justifies the use of phenomenological repulsive interactions in the multi chain 
slip-link simulations. 



P 




(41) 



of strong slip-linking (Z — > 2) the pressure approaches to zero. 



P ^ {Z ^2) 



(42) 
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Here we comment on the stress-optical rule (SOR) - in the multi chain slip-linked mod- 
els^. The rheological properties of multi chain slip-linked models are usually calculated 
by using the SOR. Namely, the stress tensor of the system is assumed to be determined 
from the average conformation tensor of subchains between slip-links (or between a slip-link 
and a chain end). In our model, these subchains clearly follow the ideal Gaussian statis- 
tics. The repulsive potential between slip-links do not alter the chain statistics. Even if 
we introduce the repulsive potential between monomers, the interaction between monomers 
are expected to be screened if the monomer density is sufficiently large and the interaction 
range is sufficiently short-"^"^. Therefore, the chain statistics is safely assumed to be Gaus- 
sian. If we ignore the contribution from the effective repulsive interaction between slip-links 
or monomers, the SOR clearly holds. However, if we consider the contribution from the 
repulsive interaction to the stress tensor, the situation is not simple. Here we consider the 
stress tensor of the system, including the contribution from the repulsive potential between 
monomers. The stress tensor of the system can be written as follows. 



-E 



^ ^^3kBT dRj{s) dR{s) 



V ^ .In b'^ ds ds 

-N 



, f J J , dUoionomeriRjis) — Rj'{s')) . , , r> / /\\ 



(43) 



where Rj represents the conformation of the j-th chain and f/monomer(^") is the repulsive 
potential between monomers (we assume that ?7monomcr('") does not diverge at r — i- 0). The 
second term in the right hand side of eq ( H3l) violates the SOR. Unfortunately, it is not 
easy to evaluate how significant this term is, because it depends on several factors such 
as the slip-linking strength or the dynamics model. Recent simulation results show that 
the rheological properties (such as the shear relaxation modulus) are not sensitive to the 
strength of repulsive interaction between monomers, except for the short time scale region—. 
This implies that the contribution from the monomer interaction to the stress tensor decays 
rapidly to its equilibrium value. Then, we expect that eq can reasonably approximate 
as the following coarse-grained form, if the slip-linking effect is not strong and the considered 
time scale is longer than the equilibration time of subchains, Tg. 

- - F E - - R,,) - ksTB, (^)' 1 (44) 

where Sj^k is the index of the A;-th slip-link on the j-th chain (we assume Sj^i < 2 < 
Sj 3 < . . . ), Rj^k = Rj{sj^k) is the position of the slip-linked point, and 1 is the unit tensor. 
Therefore, even if we introduce repulsive interactions between monomers, we can still utilize 
the SOR (at least as a good approximation). If we introduce the effective repulsive potential 
between slip-links instead of the potential between monomers, the stress tensor becomes 



-E 



^ ^^3kBTdRj{s)dRj{s) 



V ^ In 6^ ds ds 



1 dUshp-lmk{Rj,k — Rj 



(45) 

{Rj,k — Rj',k') 
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where t/siip-iink('") is the repulsive interaction between shp-hnked points. Again, the second 
term in the right hand side of eq ( H5|) violates the SOR. However, we expect that positions 
of slip-links will relax into their local equilibrium positions at the time scale longer than Tg. 
Then for such long time scales eq (H5ll will be approximated as eq (Hlj), with B2 and in 
the second term in the right hand side replaced by B2 and Z, respectively, and the SOR 
becomes approximately valid. (Of course such a naive expectation may be incorrect, and 
if so, the use of the SOR will not be justified.) From the discussions above, we consider 
that the use of the SOR in multi chain slip-link simulations is practically applicable, except 
for rather short time regions. (We do not discuss the rheological properties of our model 
explicitly, since the rheological properties depend on the dynamics model and it is beyond 
the scope of this work.) 

Before we end this section, we shortly comment on cross-linked Gaussian chain sys- 
tems^"—. We expect that our analysis will be also informative to understand statistical 
properties of cross-linked chains. In cross-linked ideal chain systems, chains are connected 
by cross-links which do not move along chains. To calculate the equilibrium statistics, the 
distribution function (or statistics) for an index along the chain, s, should be specified. By 
using the distribution function for an index, we can calculate the equilibrium statistics in the 
same way as the slip-linked chains. We expect that the equilibrium statistics of randomly 
cross-linked chains are qualitatively similar to one of slip-linked chains. Then cross-linked 
chains are expected to form compact aggregate like structures. In fact, it is known that 
cross-linked chains collapse into compact structures, without fixing some points on the sur- 
face or introducing periodic boundary conditions^"—. Our model may be utilized to study 
several Gaussian chain systems with various links which bind chains. 

CONCLUSION 

In this work we calculated the expressions of free energy and pressure for a weakly slip- 
linked multi chain system. Although these results may look physically unnatural, they are 
consistent with earlier works which pointed unusual statistics or aggregation (localization) 
behaviors in slip-linked systems. We showed that the free energy or the pressure (the equation 
of state) are different from ones for the corresponding ideal system. This is in contrast to 
the equilibrium statistics of the single chain model, which is essentially equivalent to an ideal 
chain. Unusual statistics comes from the nature of multi chain slip-link models. 

By considering the obtained expression for pressure in detail, we concluded that slip-links 
cause the effective interaction between chains. If the slip-linking effect is sufficiently strong 
chains spontaneously aggregate into compact structures. This is an artifact of the model, 
and thus we need to introduce repulsive interactions between monomers or between chains to 
avoid the unphysical aggregation. This supports previous reports about multi chain slip-link 
simulations. 

Although our analysis is mainly limited for some simplified cases, we believe that our 
results are qualitatively correct even for strongly slip-linked systems. This work will provide 
a new strategy to analyse or improve equilibrium statistics of existing multi chain slip-link 
models. 
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APPENDIX A: SINGLE CHAIN MODEL 

In this appendix, we briefly show the results for a weakly slip-linked single chain model 
We consider equilibrium statistics of slip-linked chains with a sort of mean field type approx- 
imation. In single chain models, a slip-link does not bind two chains. It is just placed on a 
polymer chain to fix it spatially. As the multi chain model case, we assume that slip-links 
are non-interacting. Then, slip-links are treated as one dimensional grand canonical ideal 
gas on a polymer chain. The effective chemical potential e is used to control the number of 
slip-links. (The physical meaning of e is different from e in the multi chain model.) If the 
slip-linking effect is sufficiently weak, we can assume that the number of slip-links on a chain 
is either or 1. Then the single chain partition function Z can be written as 

Z = ^(l + e^-/'=.T) (4g) 

The partition function of a many chain system can be calculated easily by using the single 
chain partition function (H6|) . In the foUowings, we also call the ensemble of slip-linked chains 
with the mean field type approximation as "the slip-linked single chain model" or "the single 
chain model", because its statistics is essentially determined by the single chain partition 
function (H6l) (although strictly speaking, this expression will not be appropriate). For the 
canonical ensemble of chains, the partition function of the system is expressed as 

^"M!^ "M! VAV 

where M is the number of polymer chains in the system. Thus the free energy and pressure 
can be expressed as follows 



T = knTM 



3 



In ^ - ln(l + e'/'-^) - 1 



(48) 



P = (49) 

The pressure of the single chain model is thus independent of the effective chemical potential. 
From eq ([7]), we find that the pressure of a slip-linked single chain model is the same as the 
pressure of an ideal system. 

The average number of slip-links on a chain, Z, is given by 

7 _ 1 _ 1 
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By inverting eq (150|) . the effective chemical potential e can be related to Z as 



z 



Z + 0{Z^) 



(51) 



1 - z 



where we used Z <t^l. Finally we have the following expression for the free energy. 



This expression for the free energy flS^ looks the same as the free energy of the multi 
chain model (1221) except for the numerical coefficient (1/2). However, in the single chain 
model, the average number of slip-links on a chain, Z, is independent of the chain density 
M/V. Therefore the pressure (H9l) is exactly the same as the pressure of ideal chains. This 
means that in the single chain model, the static thermodynamic properties are the same as 
the corresponding ideal system. It should be noted that the pressure is never affected by 
slip-links even if we consider strongly slip-linked cases {Z 3> 1). 

Intuitively, the differences between the single and multi chain models come from the 
difference of the slip-linking effects. Slip-links in the single chain model do not bind chains 
while slip-links in the multi chain model do. This nature is independent of the slip-linking 
strength, and thus the discussions in this appendix also holds qualitatively for strongly slip- 
linked systems. 

We should note that even in multi chain systems, the situation becomes almost the same 
if there are no direct spatial coupling between two slip-linked chains. For example, in the 
dual slip-link model by Doi and Takimoto^ or in the slip-spring model by Likhtman— , slip- 
linked two points are coupled "virtually" and not directly coupled in space. The statistics 
of a chain is not affected by slip-links. Therefore in these models slip-links do not alter 
equilibrium statistics of individual chains, and the resulting expressions of free energy or 
pressure (or other static thermodynamical properties) are essentially the same as eqs fl52|) 
and fj49|) . In a sense, we can interpret these virtual coupling type models as single chain 
models. 



We can construct clusters which contain self-slip-links (slip-links which constraint the same 
chain). As mentioned in the main text, the self-slip-linked clusters do not contribute to 
equilibrium statistics. In this appendix, we show that self-slip-linked can be neglected safely. 
We calculate the statistical probabilities for self-slip-linked and show that self-slip-linked 
spontaneously shrink and finally disappear in equilibrium. 

First Order Self-Slip-Linked Cluster 

We start from the most simple self-slip-linked structure. We consider the first order self-slip- 
link, which is a single isolated loop like structure (Figure [3]^a)). We express the indices of 
slip-hnked points by si and S2, and assume si < S2. The partition function can be formally 



jr = jr(o) _ keTKiZ + 0{Z^) 



(52) 



APPENDIX B: EFFECT OF SELF-SLIP-LINKS 
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expressed as 



ZH2) 



Q f rN rN 

dridr2 / dsi I ds2 

Jo J SI 



A6Ar2 _ 

VQ 

A3Ar2 



3A^ 



dsi / ds2 

Jsi 



3A' 



2TI{S2-Si)h\ 
-1 3/2 



3/2 



exp 



3(ri - r^)^ ^ 
■2(s2-si)62 



27r(s2 - 



(53) 



(As before, the subscript 1(2) in eq f l53|) indicates that the cluster is composed of one chain 
and two points are shp-hnked.) However, the integral in eq f l53|) does not converge, due 
to the divergence of the integrand at si = S2- Then we may interpret that the expression 
fl55]) itself has no physical meaning. Instead of the partition function flS5]) . we consider the 
probability for given si and S2, Vi(2){si, S2). This probability is proportional to the integrand 
of the last integral eq ( 153|) . 



^1(2) (si, S2) oc j dridr2 



3A' 



27r(s2 -si)62 



3/2 



exp 



3(ri - r-2)^ 
■2(s2-si)62 



A^(5(r2 - ri) 



(54) 



oc 



>2 - s,)y^ 



From eq fl54|) we obtain the following effective free energy for si and S2. 

3 

^i(2)(si,S2) = -A;sTlnPi(2)(si, S2) = -A;ijTln(s2 - si) + (const.) 



(55) 



The effective free energy (^5^ diverges to negative infinity at Si = S2, and thus we know 
that this system has no lower bound for the free energy. The system will be trapped at 
the state Si = S2 and cannot escape from that state. In other words, this system has no 
thermodynamic equilibrium state and therefore the partition function f l53|) has no physical 
meaning. Since the state Si = S2 corresponds to unentangled chain, we conclude that the 
self-slip-link is spontaneously eliminated. (These results are essentially the same as the 
theories by Sommer— , or by Metzler et al— .) 



Second Order Self-Slip-Linked Clusters with One Chain 

We proceed to more complicated cases. Here we consider the second order self-slip-linked 
structures which are composed only by one chain. There are three different second order 
self-slip-linked structures; (1) two isolated loops (Figure [3]^b)), (2) nested loops which has a 
shape of "8" (Figure [21^c)), and (3) fused loops which has a shape of "6'" (Figure El^d)). We 
may index clusters shown in Figure [3]^b), (c), and (d) as 1(4), 1(4'), and 1(4"), respectively 
(primes are introduced to distinguish clusters which have the same number of slip-linked 
points). 

For the case of two isolated loops, the situation is almost the same as the case of an 
isolated loop. There are four slip-linked points. We write indices as Si, 82,83, and S4. As 
before, we assume < Si < S2 < S3 < S4 < A^. The probability distribution for indices of 
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slip-linked points can be expressed as 

^i(4)(si,S2,S3,S4) oc J dridr2dr2,drA 



X exp 

l6 



3(r2 - ri)2 3(r3 - rs)^ 3(r4 - ra)^ 



X A^(5(r2-ri)(5(r4-r3) 
1 



(56) 



oc 



and the effective free energy becomes 

3 

^1(4) (si, S2, S3, Si) = 2^^-^ ^"^{^2 - Si) + ln(s4 - S3)] + (const.) (57) 

The effective free energy diverges at S2 = si and/or S4 = S3. This means that the system will 
be trapped at the state S2 = Si and S4 = S3, and thus the loops are spontaneously eliminated. 
The current argument can be generalized easily to many isolated loops, and thus all isolated 
loops are spontaneously eliminated. 

For the case of "8" -shaped loops, the situation is also similar. We can calculate the 
probability distribution of {sj} as 

7'i(4')(si,S2,S3,S4) oc j dridr2dr2,dri 



X exp 

,6 



[(s2-Sl)(s3-S2)(s4-S3)p/2 

3(r2 - ri)2 3(r3 - r2)^ 3(r4 - r^f 



2(S2 - Si)62 2(S3 - S2)62 

X M^5{r^ - r2)6{ri - n) 
1 



2(S4-S3)62 



(58) 



oc 



[{Sz - S2)(S4 - S3 + S2 - Si)]3/2 

The effective free energy is expressed as 

3 

^1(4') (si, 52, S3, S4) = 2^bT [ln(s3 - S2) + ln(s4 - S3 + S2 - Si)] + (const.) (59) 

As before, the effective free energy diverges at S3 = S2 and/or S4 = S3 — S2 + si. Thus the 
system will be trapped at the state S3 = S2 and S4 = si. Thus we find that the "8"-shaped 
loops are also spontaneously eliminated. 

Finally we consider the "6'" -shaped loops case. Although the method itself is the same 
as previous cases, the calculation becomes a bit complicated for this case. We have the 
following probability distribution for indices. 

Pi(4»)(si,S2,S3,S4) oc / dridr2dr3dr4 j- -r— 

J [(^2 - Si)(S3 - S2)(S4 - S3)J'^/^ 

3(r2-ri)2 3(r3-r2)2 3(r4 - ra)^ 



X exp 



2(S2 - Si)62 2(S3 - S2)b^ 2(S4 - ss)b^ 



(60) 



X A^(5(r3 - ri)5{r4 - r2) 



1 



oc 



[(S2 - Si) (S3 - S2) + (S2 - Si)(s4 - S3) + (S3 - S2)(S4 - S3)]3/2 
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and the effective free energy becomes 

3 

J'l(4//)(Si, S2, S3, S4) = -/CijTln [(S2 - Si)(s3 - S2) + (S2 - Si)(s4 - S3) + (Sg - S2)(S4 - S3)] 

+ (const.) 

(61) 

The free energy diverges at 

(S2 - Si)(s3 - S2) + (S2 - Si)(s4 - S3) + (S3 - S2)(s4 - S3) = (62) 



However, unhke the previous cases, it is not clear whether the condition ( 1621) corresponds to 
the unentangled state or not. To reduce the degrees of freedom, we try to integrate eq fl60|) 
over S2. 

dS2 7^1(4") ( Si, S2, S3, S4) 

1 (63) 

[4(S4 - S3) + (S3 - Si)]a/(s4 - S3)(S3 - Si) 

The corresponding effective free energy is 

J'i(4//)(si,S3,S4) = ^A;ijr[21n[4(s4 - S3) + (S3 - si)] + ln(s4 - S3) + ln(s3 - si)] ^^^^ 
+ (const.) 

This free energy diverges at S3 = Si and/or S4 = Si. Thus we find that the system will be 
trapped at S4 = S3 = Si. Similarly, we have the following expression for the effective free 
energy by integrating eq (I6OI) over S3. 

J'i(4//)(si,S2,S4) = ^fcBT[21n[4(s2 - Si) + (S4 - S2)] +ln(s2 - Si) + ln(s4 - S2)] ^g^^ 
+ (const.) 

Eq fl65|) gives S4 = S2 = Si as the trapped state. Therefore we find that for the case of 
"6*" -shaped loops, the system is trapped at the state S4 = S3 = S2 = Si, and the loops are 
spontaneously eliminated just like other cases. 

Second Order Self-Slip-Linked Cluster with Two Chains 

There is another second order self-slip-linked structure, which is composed by two chains 
(Figure EJ^e)). This cluster has the similar structure with clusters 1(2) or 1(4), and the 
calculation is straightforward. 

We express the indices for slip-linked points on one chain as si, S2, and S3 (si < S2 < S3), 
and the index on another chain as S4. Integrating over all variables except S2 and Si, we 
have the following probability distribution function. 

7^2(3,1) (Sl, S2) OC ^^^'3^^ (66) 
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The corresponding free energy is 

^2(3,1) (si, S2) = ^^bT [3 ln(s2 - si) - ln{N - S2)] + (const.) (67) 

As before, the system is trapped at si = S2. Thus the cluster 2(3, 1) reduces to the cluster 
2(1, 1) (Figure El^a)), which is already calculated in the main text. 

From calculations in this appendix, all the self-slip-linked clusters up to the second order 
reduce to non self-slip-linked clusters. It is natural to expect that all the higher order self- 
slip-linked clusters can also be reduced to non self- slip- linked clusters. Therefore we do not 
consider the contributions from the self-slip-linked clusters. 



APPENDIX C: CHAIN CLUSTERS APPROXIMATION 



In this appendix we derive expressions for free energy or pressure with arbitrary slip-linking 
strength, by utilizing an approximation. To make the partition function tractable, we limit 
ourselves to consider only limited set of clusters, which have one dimensionally connected, 
chain-like slip-linked structures (there are no branches or loops). We may call such clusters 
as chain clusters. The n-th order chain cluster has n chains and n — 1 slip-links. The first 
order chain cluster is a free chain, and the second and third chain clusters correspond to 
2(1, 1) (Figure [2]^a)) and 3(1, 1, 2) (Figure [2]^b)), respectively. We can calculate the partition 
function of the n-th order chain cluster n(l, 1, 2, . . . , 2) as 



n(l,l,2,...,2) 



n!A3 



(68) 



It is then straightforward to calculate the grand partition function or the grand potential for 
chain clusters. 



V 

n=l 

where we defined ^ = e^^^'^^^'^^'^ Q. The number of chains M is calculated to be 

M -- 

Thus the free energy can be expressed as follows. 



(69) 



T = kBTM 



1 1 
h -e 



kBT 



^4 

Q 



(70) 



(71) 



where ^ is related to M, V" and e via eq fiTOl) . 

The number of slip-links Z and the pressure P are then calculated to be 



1 - 



1 



(72) 
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kjrM2-l 

V 2 ^ ' 

From eq fl72l) we find that Z is a monotonically increasing function of ^, and at the hmit 
of ^ — )■ oo (this corresponds to the strong shp-hnking hmit), Z approaches to 2. We also 
find that P is a monotonicahy decreasing function of Z, and thus as the slip-linking effect 
becomes strong, the pressure decreases. From eqs ([7]) and (1731) . we have eq (HTI) . As discussed 
in the main text, the pressure approaches to zero at the strong slip-linking limit. 
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Figure 1: An image of a weakly slip-linked multi chain system. Black curves and white circles 
represent polymer chains and slip-links. Only a small fraction of chains are constrained by 
slip- links. 
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(a) (b) 




(c) (d) 




Figure 2: Schematic images for polymer chains paired by shp-hnks. Black curves and white 
circles represent polymer chains and slip-links, respectively, (a) two chains paired by a slip- 
link (2(1,1)), (b) three chains coupled by two slip-links (3(1,1,2)), and (c) and (d) two 
chains paired by two shp-links (2(1, 3) and 2(2, 2)). 
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(a) (b) (c) 




(d) (e) 




Figure 3: Schematic images for self- slip-linked polymer chains. Black curves and white circles 
represent polymer chains and slip-links, respectively, (a) a chain with an isolated loop (1(2)), 
(b) a chain with two isolated loops (1(4)), (c) a chain with nested ("8" -shaped) loops (1(4')), 
(d) a chain with fused ("6'"-shaped) loops (1(4")), and (e) paired two chains with an isolated 
loop (2(3,1)). 
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